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ABSTRACT 

We present deep K s < 21.5 (Vega) identifications, redshifts and stellar masses for most of the 
sources composing the bulk of the 24 /jm background in the GOODS/CDFS. Our identified sample 
consists of 747 Spitzer/MlPS 24 /im objects, and includes ~ 94% of all the 24 /xm sources in the 
GOODS-South field which have fluxes 5'„(24/xm) > 83 /xJy (the ~ 80% completeness limit of the 
Spitzer/GTO 24 /im catalog). 36% of our galaxies have spectroscopic redshifts (mostly at z < 1.5) 
and the remaining ones have photometric redshifts of very good quality, with a median of \dz\ = 
\Zspec Zphot |/(1 + z spec ) = 0.02. We find that MIPS 24/im galaxies span the redshift range z ~ - 4, 
and that a substantial fraction (28%) lies at high redshifts z > 1.5. We determine the existence of 
a bump in the redshift distribution at z ~ 1.9, indicating the presence of a significant population 
of galaxies with PAH emission at these redshifts. The 24 /im galaxy population ranges from sources 
with intermediate luminosities (10 10 L© < L[r < 10 11 L Q ) and low-to-intermediate assembled stellar 
masses (10 9 M Q < M < 10 n MQ) at z < 0.8, to massive (M > 10 11 M©) hyper- luminous galaxies 
(Lib, > 10 12 £q) at redshifts z ~ 2 — 3. Massive star-forming galaxies at redshifts 2 < z < 3 are 
characterized by very high star- formation rates (SFR > 5OOAf / yr) , and some of them are able to 
construct a mass of sw 10 10 M© — 10 11 M Q in a single burst lifetime (~ 0.01 Gyr - 0.1 Gyr). At lower 
redshifts z < 2, massive star- forming galaxies are also present, but appear to be building their stars 
on long timescales, either quiescently or in multiple modest burst-like episodes. At redshifts z ~ 1 — 2, 
the ability of the burst-like mode to produce entire galaxies in a single event is limited to some lower 
(M < 7 x 10 10 Mq) mass systems, and it is basically negligible at z < 1. Our results support a scenario 
where star-formation activity is differential with assembled stellar mass and redshift, and where the 
relative importance of the burst-like mode proceeds in a down-sizing way from high to low redshifts. 
Subject headings: infrared: galaxies - galaxies: evolution - galaxies: statistics 



1. INTRODUCTION 

Infrared (IR) surveys provide an unbiased way of 
studying the star-formation history of the Universe. The 
Spitzer Space Telescope (Werner et al. 2004), with its 
unprecedented sensitivity, is revolutionizing our under- 
standing of IR galaxy evolution by obtaining very deep 
IR maps. The previous space IR facilities, the Infrared 
Astronomical Satellite (IRAS) and the Infrared Space 
Observatory (ISO), set the first constraints on mid and 
far-IR galaxy evolution at z < 1 (see Genzel & Cesarsky 
2000 and Franceschini et al. 2001, for a review). Deep 
Spitzer images offer the possibility of exploring star- 
formation in the high-redshift Universe, covering the gap 
between ISO sources and the still limited quantity of 
known submillimetre sources at redshifts z ~ 2 — 3 (e.g. 
Scott et al. 2002; Borys et al. 2003; Chapman et al. 
2003). 

The study of IR sources at high redshifts is useful to 
put constraints on the early stages of star and galaxy 
formation. An important cosmological issue is the de- 
termination of what kind of galaxies host star-formation 
activity at different redshifts. At z ~ 1, a mixture of 
star-forming objects is found, from massive spiral galax- 
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ics to on-going starbursts (Franceschini et al. 2003; Bell 
et al. 2005). In addition, it is known that star- formation 
can proceed on different timescales, from bursts pro- 
duced in short time periods (typically 0.01 Gyr to 0.1 
Gyr) to sources with extended periods of quiescent ac- 
tivity. The study of the evolution of the full range of 
star-forming galaxies should allow reconstruction of the 
star-formation history of the Universe. 

Several empirical approaches based on IRAS/ISO sur- 
veys have predicted the existence of a substantial pop- 
ulation of mid-IR galaxies at high redshifts. Very re- 
cently, based on the modelling of the number counts of 
Spitzer 24 /im galaxies, Lagache et al. (2004) concluded 
that the Policyclic Aromatic Hydrocarbon (PAH) emis- 
sion characteristic of the interstellar medium (Desert, 
Boulanger & Puget 1990) should remain observable in 
galaxies up to z ~ 2.5. The main PAH spectral 
features are located at rest-frame wavelengths A r f = 
3.3,6.2,7.7,8.6,11.3,12.7,16.3 and 17 /xm. If present, 
these emission lines should enter the 24 /xm filter at red- 
shifts z ~ 6.3,2.9,2.1,1.8,1.1,0.9,0.5 and 0.4, respec- 
tively. Thus, if PAH molecules already exist in high- 
redshift galaxies, they should appear in the 24 /im pass- 
band producing the selection of a substantial population 
of sources at redshifts z > 1 (Papovich et al. 2004). 

In this work we present if s -band identifications, red- 
shifts and stellar mass estimates for most of the sources 
composing the 24 /im background down to faint fluxes in 
~ 131 arcmin 2 of the Great Observatories Origins Deep 
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Survey (GOODS)/ Chandra Deep Field South (CDFS), 
achieving 94% identification completeness for sources 
with flux S v [24 fim] > 83/iJy. The mostly complete red- 
shift identification of Spitzer 24 /im sources with magni- 
tudes R < 24 (Vega) at z < 1 has already been achieved 
by Le Floc'h et al. (2005) over an extended area of the 
CDFS. Evidence of luminous infrared galaxies up to red- 
shift z ~ 2.5 has been reported by different authors (Le 
Floc'h et al. 2004; Lonsdale et al. 2004; Papovich et al. 
2005). The evolution of the mid-IR luminosity function 
and the derived SFR densities have been explored by Le 
Floc'h et al. (2005) and Perez Gonzalez et al. (2005). 
The present work, albeit based on the study of a smaller 
area, complements the former as it exploits the unique 
quality of the GOODS datasets to obtain an almost com- 
plete identification of the sources composing the mid-IR 
background and to characterize some of their most im- 
portant properties from low to high redshifts (z ~ to 
zw3-4). 

The organization of this paper is as follows. In Section 
El we summarize the datasets and catalogs used to con- 
struct the source sample analyzed in the present study. 
In Section we present the redshift distribution of the 
24 fim galaxy population in the GOODS/CDFS. In Sec- 
tion 01 we give estimates of their characteristic stellar 
masses. In Section we study the evolution of IR lu- 
minosities and stellar masses with redshift and put con- 
straints on the star-formation rates (SFR) and timescales 
for star-formation activity in different-mass galaxies. Fi- 
nally, in Sectional we summarize our results and present 
some concluding remarks. We adopt throughout a cos- 
mology with H = 70kms~ 1 Mpc _1 , J7m = 0.3 and 
fl A = 0.7. 

2. iCs-BAND COUNTERPARTS TO MIPS 24 /im SOURCES 
2.1. The data samples and redshift estimates 

Deep 24 fim observations of ~ 2000 arcmin 2 of the 
CDFS have been carried out with the Multiband Imag- 
ing Photometer for Spitzer (MIPS; Rickc et al. 2004), 
as part of the Guaranteed Time Observers (GTO) pro- 
gram. The data reduction, extraction of sources and 
photometry measurements are explained in detail in Pa- 
povich et al. (2004). The source catalog achieves 80% 
completeness at a flux limit of S v ~ 83 /iJy and the 
fraction of spurious sources at that limit is determined 
to be <10%. The 24 /mi source density in the CDFS 
is N(S V > 83/iJy)= 4.5arcmin~ 2 . Above a flux of 
S v ~ 60 fiJy, resolved Spitzer sources account for ~ 70% 
of the 24 /xm cosmic infrared background (Papovich et al. 
2004; Lagache et al. 200 4). 

The GOODS project l|Giavalisco et al. 20041) has pro- 
vided deep multiwavelength data for ~ 160 arcmin of 
the CDFS. The data products are released to the astro- 
nomical community in a fully reduced mode. As part 
of the European Southern Observatory GOODS Imag- 
ing Survey program (ESO GOODS/EIS; Vandame et al., 
in preparation), deep J and K s -h&n& images taken with 
the Infrared Spectrometer and Array Camera (ISAAC) 
on the Very Large Telescope (VLT) have been made pub- 
lic for ~ 131 arcmin 2 of the GOODS/CDFS. This is the 
area analysed in the present work. The ESO/EIS imag- 
ing in this field is complemented by deep imaging data 
in the B, V, 2775 and zsso bands taken with the Ad- 



vanced Camera for Surveys (ACS) on board the Hub- 
ble Space Telescope (HST) and by 3.6 fim, 4.5 /xm, 5.8 /im 
and 8.0 / im data taken wit h the Infrared Array Camera 
(IRAC) l|Fazio et al. 2004|) on board Spitzer. 

A catalog of K s < 21.5 (Vega) sources has been se- 
lected in the GOODS/EIS CDFS. Photometric redshifts 
based on up to eight broad-bands {BVItjsZ&zqJKs and 
IRAC 3.6 /im and 4.5 /im bands) have been obtained for 
the whole sample, using the publicly available multiwave- 
length data and the public code HYPERZ (Bolzonella, 
Miralles & Pello, 2000) with the GISSEL98 spectral en- 
ergy distribution (SED) template library (Bruzual & 
Chariot, 1993). Dust-corrections have been taken into 
account through the convolution of the SED templates 
with the Calzetti et al. (2000) reddening law. 

Although U-band data exist for the CDFS, they have 
not been included in the input catalogs for the photo- 
metric redshift algorithms as the U-band images have 
shallower depth and poorer resolution than the ACS and 
ISAAC images in the same field. However, U-band im- 
ages have been used to control the HYPERZ output: for 
every i^ s -band source with a counterpart in the shallow 
U-band catalogs, the resulting photometric redshift was 
determined as the best-fit value constrained to a max- 
imum redshift z p ) lot — 2, as higher redshift sources are 
unlikely to be bright at such short wavelengths. In addi- 
tion, the public code BPZ (Bcm'tcz 2000) has been used 
to obtain a second, independent, set of redshift estimates 
for all the K s < 21.5 sources. In the cases of sources 
with HYPERZ photometric redshifts z v hot > 2 not con- 
firmed by BPZ, we adopted the lower estimates from the 
BPZ code. In addition, pub lic COMB017 photomet- 
ric redshifts fWolf et al. 2004) have been used to replace 
the redshift estimates of those sources with magnitude 
R < 23.5 (Vega) and z < 1 (the most accurate regime 
for COMBO 17 redshifts). Further details on the redshift 
estimations for the total K s < 21.5 sample are given 
in a separate paper (Caputi et al. 2005b). Figure^,) 
compares the obtained photometric redshifts with spec- 
troscopic redshifts for all those K s < 21.5 sources in the 
GOODS/EIS CDFS included in different publicly avail- 
able spectroscopy samples (e.g. Vanzella et al. 2005; 
Le Fevre et al. 2004). We observe a very good agree- 
ment between the photometric estimates and the real 
redshifts in most cases, with a median for the absolute 
relative errors of \dz\ = \z spec - z phot \/ (1 + z spec ) = 0.02. 
Figured) shows the corresponding histogram of relative 
errors. The rms of the distribution is 0.03. However, in 
spite of the good accuracy obtained for the photomet- 
ric redshifts, the incorporation of available spectroscopic 
redshifts is of much benefit for any study of a galaxy red- 
shift distribution. Thus, in order to maximize the qual- 
ity of the final redshift catalog, the existing spectroscopic 
redshifts in the CDFS have been used to replace the pho- 
tometric estimates of the K s < 21.5 sources whenever 
possible (i.e. 23% of the total K s < 21.5 sample). 

2.2. Cross- correlation of the catalogs 

We cross-correlated the K s < 21.5 GOODS/CDFS 
source catalog with the Spitzer/MIPS 24 /im 
GTO/CDFS catalog both to investigate and put 
constraints on the mid-IR emission of ^-selected galax- 
ies and to identify and characterize the optical/near-IR 
properties of the Spitzer/MIPS 24 /im sources in the 
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GOODS South field. This work presents the results of 
the latter, while a full analysis of 24 /tm-detected and 
non-detected if s -selected galaxies in the GOODS/CDFS 
will be presented elsewhere (Caputi et al. 2005c, in 
preparation). In addition, the IRAC and MIPS proper- 
ties of galaxies selected with colours J — K s > 2.3 are 
studied by Papovich et al. (2005). 

We looked for 24 //m counterparts of the K s < 21.5- 
selected sources in the GOODS/CDFS using a match- 
ing radius of 2", to minimize the number of multi- 
ple identifications. The astrometric accuracy of the 
MIPS-24/tm images (e.g. Le Floc'h et al. 2005) allows 
the use of such a small matching distance. We found 
812 K s < 21.5-selected sources associated with MIPS 
24 /mi sources within this radius. In 65/812 cases, we 
found two or more A^-selected sources associated with 
a same 24 /im source. In these cases we selected the 
closest counterpart, leaving 747 24 /tm sources with a 
K s < 21.5 identification. We identified a negligible frac- 
tion (11/747) of galactic stars among these MIPS 24 /mi 
sources, and we exclude them from all the following anal- 
ysis. We determined redshifts for each of the remaining 
736 24 /im sources. 36% of them have spectroscopic red- 
shifts and 21% more have COMBO 17 photometric red- 
shifts. The remaining 43% have HYPERZ/BPZ photo- 
metric redshifts. 521/747 sources are above the 24 /mi- 
catalog 80% completeness limit, i.e. have S v > 83 /tJy. 
On the other hand, we find 4 fC s -selected objects that 
could be associated with two or more 24 /im sources. 
All of these cases correspond to low redshift galaxies, 
with spectroscopic redshifts z spec < 0.13. Individual 
inspection of both the K s and 24 /jm images suggests 
the Spitzer/MIPS detection of multiple nodes of mid-IR 
emission for each of these nearby galaxies. In Figure 
we show the total 24 /im flux S v versus K s magnitude 
for the 747 MIPS sources with K s < 21.5 identifications 
in the GOODS/CDFS. In the 4 cases of multiple MIPS 
sources associated to a single if s -band object, we con- 
sidered that the 24 /tm flux was the sum of the differ- 
ent 24 /tm components. 686/747 MIPS 24 /im sources in 
the GOODS/CDFS have counterparts classified as nor- 
mal galaxies (circles in Figure|2J), while 50/747 are active 
ones (AGNs or QSOs; squares in Figure EJ. The iden- 
tification of K s sources with active galaxies is based on 
the cross-corre lation with the X-r ay catalogs available 
for the CDFS ijSzokolv et al. 2004|) . This approach pro- 
vides only a lower limit to the fraction of active galaxies 
among the 24 /tm galaxy population (Alonso-Herrero et 
al. 2005; Donley et al. 2005). The cross-correlations of 
the K s < 21.5 GOODS/CDFS source catalog with the 
Spitzer/MIPS 24 /tm GTO/CDFS catalog allows us to 
achieve a completeness limit of ~ 94% for the identifica- 
tion of MIPS 24 /im sources with flux S v > 83 /iJy in the 
GOODS/CDFS (cf. Tabled). 

As we mentioned above, our HYPERZ photometric 
redshift determinations are based on the SED fitting 
of the galaxies made using the templates in the GIS- 
SEL98 library of Bruzual & Chariot. Also, this SED 
fitting procedure has been used to determine derived pa- 
rameters (e.g. estimated stellar mass) for all the galax- 
ies, independently of whether the redshifts were from 
HYPERZ/BPZ, COMB017 or spectroscopic. The GIS- 
SEL98 library is composed of a wide range of synthetic 
SEDs based on stellar spectra but does not have any 



optical to near-IR power-law SED. To investigate the 
impact of this limitation in our redshift distribution of 
MIPS 24 /tm galaxies, we identified the presence of 24 /tm 
power-law SED galaxies in the GOODS/CDFS using the 
catalog of IR-power-law sources constructed by Alonso- 
Herrero et al. (2005). We found 20 24 /tm galaxies with 
IR-power-law spectra in the GOODS/CDFS. 9 of these 
power-law galaxies have spectroscopic redshifts, leaving 
only 11 power-law sources with a HYPERZ/BPZ photo- 
metric redshift (and 7/11 are X-ray detected; see Section 
0J). Consequently, plausible erroneous redshift estimates 
due to inadequate SED template fitting are reduced to 
a very minor fraction of the 24 /tm sources studied here 
and, thus, they should have a basically negligible impact 
on all the analysis in this work. 

2.3. Summary of properties of the 24 /tm source sample 
analyzed in this work 

For clarity, we summarize the properties of the source 
sample analyzed in this work: 

• our sample is composed of the 747 24 /im sources in 
- 131 arcmin 2 of the GOODS/CDFS which have a 
K s < 21.5 (Vega) counterpart. 521/747 sources 
have a flux above the GTO catalog 80% com- 
pleteness limit, i.e. Sy(24/im) > 83/iJy. These 
521 objects constitute ~ 94% of all the MIPS 
24 /tm sources with S v (24 /tm) > 83 /i Jy in the 
GOODS/CDFS. 

• Only 11/747 sources in our sample are identified as 
galactic stars. 50/747 sources are X-ray classified 
AGNs or QSOs. The remaining 686/747 sources 
are normal galaxies (i.e. no X-ray AGN or QSO). 

• We determined redshifts for each of the 736 galax- 
ies in our sample (i.e. all sources except the 11 
galactic stars). 36% of these redshifts are spectro- 
scopic. 21% additional redshifts have been taken 
from the COMB017 photometric redshift catalog 
(only for sources with R-band magnitude R < 23.5, 
Vega, and redshifts z < 1). The photometric red- 
shifts for the remaining 43% of the sample have 
been obtained using HYPERZ/BPZ and are based 
on broad-band photometry covering from the opti- 
cal B-band through the 4.5 /tm band. 

In this way, our sample allows us to identify and charac- 
terize the evolution down to faint fluxes for most of the 
sources comprising the 24 /tm background. 

Perez-Gonzalez et al. (2005) have carried out a study 
of MIPS 24 /tm detections using a different photometric 
redshift technique. Their approach uses empirical SEDs 
rather than synthetic templates and it can fit redshifts 
using a wider range of multiwavelength data. Their work 
has the advantage of having been applied to larger fields. 
However, high signal-to-noise data may not be available 
in all bands for these large areas and thus, their resulting 
redshifts have a typical accuracy of about 10%, signifi- 
cantly worse than those obtained here. Our approach de- 
pends on high-quality and homogeneous input data, only 
available for very limited areas of the sky. By combining 
the results of both works it is possible to: 1) determine 
overall trends and the effects of cosmic variance based 
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on the approach of Perez- Gonzalez et al.(2005); and 2) 
prove the redshift distribution more accurately in a more 
limited sky region based on the approach in this paper. 

3. THE REDSHIFT DISTRIBUTION OF GALAXIES IN 

THE GOODS/CDFS 

Figure^) shows the redshift distribution of the MIPS 
24 /im galaxies with K s < 21.5 counterparts in the 
GOODS/CDFS. The empty (shaded) histogram corre- 
sponds to all (S v > 83/iJy) sources. The peaks in the 
redshift distribution at z ~ 0.7 and z ~ 1.1 are due to 
the effect of large-scale structure in the CDFS (e.g. Le 
Fevre et al. 2004). We find that 72% of the 24 /im galax- 
ies lie at redshifts z < 1.5, while the remaining 28% are 
found to be at z > 1.5 (the percentages are similar for the 
samples containing all and only S v > 83 /iJy galaxies) . 
This confirms the existence of a substantial population of 
mid-IR sources at high redshifts. Approximately half of 
the z > 1.5 sources have the characteristic colours of Ex- 
tremely Red Galaxies (ERGs), (J 775 -K s ) > 4.0 (Vega), 
indicating that an important fraction of the high-redshift 
mid-IR background is constituted by optically obscured 
sources, as possibly expected (cf. Yan et al. 2004). 

Figure (3t>) shows the normalized redshift distributions 
of the MIPS 24 /im galaxies with K s < 21.5 counterparts 
(black dashed and solid lines for all and S v > 83 /iJy 
sources, respectively), compared to the redshift distribu- 
tion of the total K s < 21.5 galaxy population (grey solid 
line) in the same field. Several features are present in 
all the three curves, which are the consequence of cosmic 
variance effects. In contrast, we observe the presence of a 
depression in the redshift distribution of 24 /im galaxies 
at redshift z ~ 1.5 and a bump at redshift z ~ 1.9, both 
of which do not appear in the total K s < 21.5 galaxy 
curve. 

To assess the significance of the features observed 
exclusively on the 24 /im galaxy redshift distributions, 
we computed confidence limits on the different curves. 
We performed Monte-Carlo simulations to create 1000 
mocked redshift catalogs, alternatively for the total K s < 
21.5 galaxy population and for the MIPS 24 /im galaxies 
with flux S v > 83 /iJy. We constructed each mock cat- 
alog assigning to each source a random redshift, with a 
gaussian probability centred at the original redshift z of 
the source and a dispersion equal to 0.02 (1 + z) (i.e. the 
median error at the corresponding redshift). The red- 
shifts of those sources with spectroscopic values were left 
fixed. We re-computed the normalized redshift distribu- 
tion for each mock catalog and determined confidence 
limits on the original distribution curves. The 95% con- 
fidence limits on the total K s < 21.5 galaxy and the 
MIPS (S v > 83 /iJy) galaxy distributions are shown in 
Figure 0] 

Figure 0] shows that, even taking into account the er- 
ror bars, the peak in the 24 /im redshift distribution at 
z ~ 1.9 is significant, lying ~ 4cr from the original to- 
tal K s < 21.5 galaxy curve. Thus, we conclude that 
the redshift distribution of 24 /im galaxies presents a 
real secondary bump at these high redshifts. This sec- 
ondary peak has been predicted by Lagache et al. (2004) 
and is the consequence of the selection effect produced 
by the presence of PAH emission features entering the 
observed 24 /im band. Given the width of the 24 /im 
filter (whose transmission covers the wavelength range 



~ 20 — 28 /im), both the 7.7/xm and the 8.6/xm PAH lines 
could contribute to the redshift distribution peak ob- 
served at z ~ 1.9. Our results allow us to conclude that 
PAH molecules must be already present in star-forming 
galaxies at high redshifts. 

The confidence limits shown in Figure 01 also indicate 
that the depression observed in the 24 /im redshift dis- 
tribution at z ~ 1.5 is only marginally significant within 
our sample. The existence of such a depression could be 
interpreted as due to the 9.8 /im silicate absorption fea- 
ture entering the 24 /im filter at this redshift, which has 
been observed for some galaxies in previous works (Houck 
et al. 2005; Yan et al. 2005). However, given the errors 
in the redshift distributions presented here, we cannot 
reach any firm conclusion on the possible selection effect 
produced by silicate absorption at high redshifts. 

We also see a substantial deficit of 24 /an sources at 
redshift z ~ 0.9, with respect to the total K s < 21.5 
sources. At this redshift, a positive selection effect on 
24 /im galaxies would be expected due to the presence of 
the 12.7 /xm PAH emission line. The observed relative 
deficit of bright mid-IR selected galaxies indicates that, 
within our sample, star-forming galaxies are present in a 
minor proportion at this redshift. This is quite likely due 
to a mere cosmic variance effect. However, we note that 
this deficit in our 24 /im sample occurs at the same red- 
shift where there is a dip in the total K s < 21.5 redshift 
distribution. If we consider that galaxy interactions are a 
triggering mechanism for star-formation, then it will be 
possible that star-formation activity is particularly in- 
hibited in regions with an under-density of objects. The 
study of similar regions in other areas of the sky are nec- 
essary to determine whether the two facts are actually 
related. 

Evidence of the existence of PAH emission in the spec- 
tra of a few high-redshift 1.7 — 1.8 < z < 2.6 — 2.8 galax- 
ies has recently been presented by Houck et al. (2005), 
Yan et al. (2005) and Lutz et al. (2005). In this work, 
we extend this evidence through the study of the red- 
shift distribution of the whole 24 /im galaxy population. 
It is interesting to note, however, that both Houck et 
al. (2005) and Yan et al. (2005) found that only a minor- 
ity of their galaxies at z ~ 1.7 — 2.8 showed clear PAH 
emission features in their spectra, while most of their re- 
maining high-redshift galaxies were AGN-dominated or 
had silicate absorption features. The difference in the 
nature of the sources dominating Houck et al. and Yan 
et al. samples can be explained taking into account that 
their objects are on average an order of magnitude more 
luminous than most of the sources analyzed in this work. 
Deep 24 /im samples are necessary to uncover a substan- 
tial population of PAH galaxies at redshift z ~ 2. Fig- 
ure [5] shows the 24 /im flux of each galaxy S v versus red- 
shift z. The symbols are the same as in Figure Filled 
symbols correspond to sources with spectroscopic red- 
shifts. Within our sample, we see that only a few bright 
S v (24fim) > 0.5 mJy sources are placed at high redshift 
(z > 1.5) and they are all active galaxies. At fainter 
fluxes S v < 0.3 mJy, starbursts produce the bulk of the 
mid-IR emission at high redshifts. 

Perez-Gonzalez et al. (2005) found that 24% of the 
24 /im galaxies with S u > 83 /iJy in an extended region 
of the CDFS and the Hubble Deep Field North (HDFN) 
were at redshifts z > 1.5, in agreement with our value of 
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28% within the errors and taking into account the differ- 
ent fields surveyed. However, their redshift distribution 
does not show the PAH-induced peak at 1.6 < z < 2.2. 
They suggested that this feature may have been blurred 
by the errors in photometric redshifts; this suggestion is 
confirmed by our detection of the PAH bump using red- 
shifts of higher accuracy. Taken together, the two stud- 
ies show that the drop in the number of sources beyond 
z ~ 1.2 holds generally and is not just observed in the 
CDFS. In addition, we show that there is a significant 
peak produced by PAH emission on top of this general 
trend. 

Figure shows the comparison of our observed red- 
shift distribution of 24 /jm galaxies in the GOODS/CDFS 
(solid histogram) with the distribution predicted by La- 
gache et al.(2004) using the Lagache, Dole & Pugct 
(2003) model (dashed histogram). Both distributions 
correspond to sources with flux S v > 83 /zJy in an area 
of 131 arcmin 2 . The observed and predicted distribu- 
tion have 521 and 518 galaxies, respectively, confirming 
that the Lagache et al. (2003) model very well repro- 
duce the observed number counts. Also in agreement 
with the model, we find that a significant population of 
24 (im galaxies lie at redshifts z > 1.5, and our accurate 
redshift determinations reveal the predicted secondary 
peak in the redshift distribution at z ~ 1.9 produced 
by the presence of PAH emission. However, the fraction 
of galaxies at redshifts z > 1.5 is significantly overpre- 
dicted by the Lagache et al. model. Other models (e.g. 
Chary & Elbaz 2001) basically predict a negligible frac- 
tion of sources at redshifts z > 1.5 and also do not fit 
the observed redshift distribution. Therefore, theoreti- 
cal models require further refinements to reproduce the 
distant infrared-emitting galaxy populations. 

4. THE ASSEMBLED STELLAR MASSES OF MIPS 24/zm 
GALAXIES 

Figured shows the estimated assembled stellar masses 
of the MIPS 24 /j,m galaxies versus redshift z in the 
GOODS/CDFS. The estimation of the stellar masses is 
based on the optical-to-near-IR SED fit of each galaxy 
at the determined redshift (either spectroscopic, from 
COMB017 or HYPERZ/BPZ) and is completely inde- 
pendent of its 24 /mi properties. We determined esti- 
mated stellar masses for all the 24 /xm galaxies with a 
K s < 21.5 identification, except for any X-ray-classified 
AGN/QSO or any other X-ray source without a suit- 
able GISSEL98 template optical-to-near-IR SED fit. We 
excluded the latter objects to avoid galaxies with par- 
tial contamination of the optical-to-near-IR light by a 
hidden AGN. The stellar masses are computed from the 
modelled rest-frame i"T s -band galaxy luminosity of each 
galaxy, where the mass-to-light ratios have a minimum 
dependence on the SED star-formation histories or dust 
corrections (see Caputi et al. 2005a, b for further details). 
The resulting stellar mass estimates are typically accu- 
rate within a factor < 2. A single power-law Salpeter 
IMF over stellar masses M = (0.1 — 100) M Q has been 
assumed throughout. The estimated mass-completeness 
limits for counterparts of the MIPS 24 /im galaxies are 
1.5xl0 10 M Q , 7.0xl0 10 M Q andl.5xlO n M at redshifts 
z = 1,2,3, respectively, based on the K s — 21.5 magni- 
tude limit and the median of the k-corrections. However, 
it should be noted that the mass completeness limits are 



basically irrelevant for the sources with flux S v > 83 /iJy 
(symbols filled with a cross in Figure 0), as the i^-band 
identifications are almost complete above this limit. 

Inspection of figure shows that, at least at redshifts 
z < 1, star-formation activity takes place in galaxies of a 
wide range of assembled stellar masses, from ~ 10 9 M Q 
to ~ 1O 12 M . At high redshifts z > 2.5, the existence 
of a significant population of massive M > 10 11 M Q star- 
forming galaxies is revealed in the mid-IR at the depth 
of the GTO/CDFS images. It is now commonly believed 
that a fraction of present-day massive galaxies is already 
in place at high redshifts (Cimatti et al. 2004; Glaze- 
brook et al. 2004; Caputi et al. 2005a, b). However, it is 
still unclear when the bulk of the stars in these massive 
galaxies was formed (cf. Papovich et al. 2005). The iden- 
tification of massive objects on the MIPS 24/im images 
can put constraints on their star-formation histories and 
the amount of stellar mass built-up at different redshifts. 

5. CONSTRAINTS ON THE STAR-FORMATION HISTORY 
5.1. The evolution of star-forming galaxies with redshift 

The study of the physical properties of IR-selected 
galaxies at different redshifts is of importance to under- 
stand how star-formation evolved with cosmic time. In 
this section, we use the 24 fluxes of the galaxies in 
the GOODS/CDFS in conjunction with the redshift and 
stellar mass estimations, to study the evolution of star- 
formation activity. 

Chary & Elbaz (2001) and Elbaz et al. (2002) showed 
that the mid-IR luminosities of nearby galaxies were cor- 
related (with some scatter) with their bolometric IR lu- 
minosity Lju, defined as Ljr — L(8 — 1000/im). They 
fitted the following set of relations: 

L IR = ll.l+_l 5 7 x (^[15/im]) - 998 , (1) 

= 0.89t°;p x (uL4l2^m]) lm \ (2) 

= 4.78± 2 ;g x (^[6.75/H)°- 998 (3) 
for uL u [6.75^m] < 5 x lO 9 ^, 

= 4 - 37 -2'.i3 x 10_6 x (*aM6-75aH) 1,62 (4) 
for vL v [6.75nm] > 5 x 10 9 L Q , 

where all the luminosities are in solar units. There are 
many indications that the overall spectral energy distri- 
butions of high-redshift infrared galaxies are similar to 
the local ones captured in equations (1) to (4). For ex- 
ample, the local far-IR/radio and mid-IR/radio correla- 
tions (Condon 1992) still hold at higher redshifts (Ap- 
pleton et al. 2004). Composite SEDs of high redshift 
infrared galaxies resemble local templates (e.g. Egami et 
al. 2004). The MIPS observations in the CDFS indicate 
that, at least to z ~ 1, galaxies have 24 and 70 /im flux 
ratios that follow the local distribution as a function of 
total IR luminosity. In addition, number-count models 
assuming similar behaviour can fit the 24, 70, and 160/im 
data simultaneously (Lagache et al. 2004). 

We used the 24 ^m fluxes of the MIPS galaxies with a 
K s < 21.5 counterpart in the GOODS/CDFS to compute 
their bolometric luminosities Ljr, assuming the Chary 
& Elbaz relations in eqs. (1) to (4). We excluded all 
the known AGNs/QSOs from the present analysis. We 
considered that the 24 /zm flux mapped the rest-frame 
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15/im, 12/zm and 6.75/im fluxes in the redshift ranges 
0.4 < z < 0.8, 0.8 < z < 1.2 and 2.0 < z < 3.0, respec- 
tively. In the redshift range 1.2 < z < 2.0, we assumed 
that the infrared luminosity Lm was given by the av- 
erage of the luminosities obtained using eqs. (2) and 
(3-4). Figure |H1 shows the bolometric IR luminosities of 
the 24 /jm galaxies versus assembled stellar mass in dif- 
ferent redshift bins. The dotted lines delimit the region 
of luminosity completeness at the mean redshift of each 
bin, taking into account the «S'„(24/im) = 83 //Jy limit 
and the fact that our .ft" s -band identifications are almost 
complete above that flux (i.e. the mass completeness lim- 
its imposed by the K s = 21.5 cut are basically irrelevant 
above that limit). 

Although for clarity the error bars are not shown in 
Figure |H1 we estimated the errors on our computed bolo- 
metric luminosities using an independent set of models. 
We integrated the different Dale et al. (2001) and Dale & 
Helou (2002) SED templates, normalized to our observed 
24 /xm fluxes, to obtain the corresponding IR luminosi- 
ties as a function of redshift. The complete set of the 
Dale et al. SEDs consists of 64 templates, each one char- 
acterized by a parameter a (0.0625 < a < 4.0), which 
depends on the combination of the observed local SEDs 
used to construct these models. Compared to the values 
obtained with the Dale et al. templates, our luminosi- 
ties based on the Chary & Elbaz formulae appear to be 
accurate within a factor 2-3 to redshift z « 1.5. Above 
that redshift, the consideration of the complete set of 
the Dale et al. SEDs would yield uncertainties of up 
to a factor 10-20 (Dale et al. 2005). However, Dale et 
al. (2005) recently used the Dale et al. templates to 
fit the observed 1-850 fim photometry of a sample of 75 
local galaxies. The majority of these galaxies ( > 75%) 
were best-fitted by templates with 2.0 < a < 4.0. As 
we explained above, there are many indications that the 
overall spectral energy distributions of high-redshift IR 
galaxies are similar to the local ones. Thus, this subset of 
the Dale et al. templates should be representative of the 
majority of the IR galaxy types also at high redshifts. 
If we restrict the comparison of our luminosities com- 
puted with the Chary & Elbaz formulae to those values 
obtained with the Dale et al. 2.0 < a < 4.0 templates, 
we find that the error bars would still be within a factor 
2-3 up to redshift z rs 2, and a factor 5 would be suffi- 
cient to account for the errors above that redshift. These 
factors should be representative of the error bars on the 
IR luminosities of the majority of our galaxies, including 
the errors introduced by the k-corrections within each of 
the considered redshift bins. 

We observe that the evolution of the IR luminosities 
with redshift is still very significant even taking into ac- 
count the error bars. Panel a) shows that most of the 
24 /im galaxies at redshifts 0.4 < z < 0.8 have infrared 
luminosities Lir < 10 11 Lq. The maximum observed in- 
frared luminosities increase with redshift, and luminous 
infrared galaxies (LIR Gs) characterized by 10 11 L Q < 
L IR < 10 12 L Q USanders fc Mirabel 1996) are the dom- 
inant 24 /zm sources at redshifts 0.8 < z < 1.2 at the 
depth of the CDFS MIPS/GTO images (Le Floc'h et al. 
2005). The majority of the mid-IR sources at 0.4 < z < 
1.2 are hosted by intermediate-mass galaxies with stel- 
lar masses 10 10 Mq < M < 10 11 Mq, in agreement with 



recent findings by Hammer et al.(2005), although some 
more massive galaxies could also be classified as LIRGs 
at these redshifts. Within our surveyed area, there is vir- 
tually no ultra-luminous infrared galaxy (ULIRG) with 
Lin > 10 12 Lq at z < 1.2. ULIRGs might be present at 
these low redshifts, but are indeed very rare (e.g. Flores 
et al. 1999). At 1.2 < z < 2.0, ULIRGs start to be a very 
significant population (~ 65% at S v > 83/iJy) and they 
are also hosted by intermediate to high-mass galaxies. At 
z > 2, the limits of the 24 /jm survey only allow us to ex- 
plore the bright IR-luminosity end of the star-formation 
activity. We observe sources with extremely high infrared 
luminosities 10 12 Lq < Lju < 10 14 Lq mainly harboured 
by galaxies with stellar masses M > 10 11 Mq. Our re- 
sults on luminosity and SFR evolution agree well with the 
studies for larger fields made by Le Floc'h et al. (2005) 
and Perez-Gonzalez et al. (2005). 

Thus, figure |H1 confirms the global trend of a high de- 
gree of evolution in the star-formation activity from low 
to high redshifts. The extremely high IR luminosities 
observed at z > 2 show that star-formation has been a 
much more violent process in the past and these violent 
episodes of star- format ion must have built up a signifi- 
cant part of the stellar mass content of at least a fraction 
of present-day massive galaxies. We explore in more de- 
tail this possibility in the next section, where we analyze 
and constrain the typical timescales for star-formation of 
galaxies with different assembled stellar masses at differ- 
ent redshifts. 

5.2. Probing star-formation in different-mass galaxies 

An estimator of the instantaneous star-formation rate 
(SFR) of a galaxy based on its bolometric infrared lumi- 
nosity has been obtained by Kennicutt (1998): 

SFR = 1.72 x lO" 10 ^, (5) 

where the star- format ion rate SFR is given in Mq/yt, 
the luminosity Lju is in solar units and a Salpeter IMF 
over stellar masses M = (0.1 — 100) Mq is assumed. We 
used equation (JSJ) to obtain an estimate of the instan- 
taneous SFR in the 24 /im galaxies studied here, and to 
explore the relation of this SFR with the already assem- 
bled stellar mass of each galaxy, at different redshifts. 

Figure shows the evolution of the SFR with redshift 
for galaxies with assembled stellar masses in different 
mass bins. In all the panels, the dashed line delimits the 
region of completeness of the observed SFR versus red- 
shift, as imposed by the flux limit S^(2Afim) = 83/iJy. 
We find that the IR-derived SFR exponentially grows 
with redshift, independently of the galaxy stellar mass 
(of course the lower envelope of SFR is a consequence 
of the limits of the 24 /im survey, but the maximum ob- 
served values do exponentially grow with redshift). How- 
ever, panel a) shows that, at a given redshift, most of the 
lowest mass galaxies have relatively small SFR, which 
are at most marginally above the completeness limit of 
the 24 /jm survey analyzed here. This effect is particu- 
larly evident at high redshifts z > 2, where we observe 
a lack of galaxies with stellar mass M < 1.5 x 1O 1O M0 
and SFR above the completeness limit. Unless the minor 
fraction of unidentified sources with S , „(24/im) > 83/iJy 
constitutes a biased population of low mass galaxies at 
high redshifts, we conclude that low-mass objects with 
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very high SFR (SFR > 5OOM /yr) are indeed rare at 
z > 2. On the contrary, galaxies with greater stellar mass 
M > 1.5 x 10 10 Af Q and SFR > 500M Q /yr are progres- 
sively found at higher redshifts. 

Figure EH shows the ratio of already assembled stel- 
lar mass over the instantaneous SFR, Al/SFR, versus 
redshift z for the 24 (im galaxies with a K s < 21.5 coun- 
terpart in the GOODS/CDFS. Different symbols corre- 
spond to galaxies with stellar masses in different mass 
ranges, as in Figure The solid line shows the age of 
the Universe as a function of redshift, while the dashed 
line indicates the upper limit of the typical lifetime of 
a starburst (~ 0.1 Gyr). The assembled stellar masses 
have been computed as in Section^ and are completely 
independent of the 24 fim properties of each galaxy. If the 
SFR in a galaxy were larger in the past or at most stayed 
constant over cosmic time, the ratio Mj SFR would give 
an upper limit for the lifetime of the source. Therefore, 
figure lTUI allows putting constraints on the star-formation 
history of the MIPS 24 /im galaxies present at different 
redshifts. 

A first striking conclusion is that, at high redshifts 
z ~ 2 — 3, a starburst lifetime appears to be suffi- 
cient to construct a substantial fraction and up to the 
whole amount of the stellar mass of some massive galax- 
ies. If equations (1) to (5) are still applicable at those 
redshifts, very high SFR are predicted, ranging from 
~ 500A/ Q /yr to 5000M©/yr. In this way, a stellar 
mass of ~ 10 10 Mq — 10 11 M Q could be constructed in 
a ~ 0.01 Gyr — 0.1 Gyr period of time. However, this 
does not necessarily mean that all these systems are be- 
ing formed at z ~ 2 — 3. Recent works suggest that 
a non-negligible fraction of present-day massive galaxies 
would be present before even higher (z > 3) redshifts 
(van Dokkum et al. 2004; Caputi et al. 2004; Caputi 
et al. 2005a, b). Certainly, we find other 24 /im mas- 
sive galaxies for which the derived lifetimes Mj SFR are 
larger (up to ~ 1 Gyr), indicating that these galaxies 
would be present since higher redshifts and that star- 
formation histories longer than a single burst lifetime 
are necessary to build the stellar mass present in them. 
Moreover, other massive galaxies with SFR below the 
limits of the 24 /im survey analyzed here also exist at 
redshift z ~ 2 — 3 (Caputi et al. 2005c, in preparation). 

Interestingly, we note that the best-fit optical to near- 
IR SEDs for our 24 /im galaxies at redshifts z > 2 cor- 
respond to different galaxy ages, many of which exceed 
the maximum ~ 1 Gyr lifetime derived with the ratio 
Mj SFR. This fact could be indicating that the instanta- 
neous star-formation rates derived from the 24 /im fluxes 
of these galaxies were not constant or larger in the past. 
In order for these galaxies to have ages of ~ 1-3 Gyr 
(as those derived from the best-fit SEDs), their star- 
formation rates must have been smaller at some moment 
in the past. This is strongly suggesting that some mas- 
sive galaxies could be already in place at higher redshifts 
and experience a non-continuous star-formation history, 
with temporary episodes of high star-formation activity 
(cf. also Papovich et al. 2005b). Those suffering such 
episodes at z ~ 2—3 would be detected as hyper luminous 
infrared galaxies at these redshifts in mid-IR surveys. If 
this is the case, the major episodes of star-formation at 
z ~ 2 — 3 would not account for the complete assembly 



history of these massive galaxies. A similar argument 
has been given by Hammer et al. (2005), who propose 
that the stellar mass growth of intermediate-mass galax- 
ies at z < 1 could proceed by successive star-formation 
episodes characterized by a LIRG phase. 

As we mentioned in Section|3 a substantial fraction of 
the 24 /xm sources at high redshifts have the characteris- 
tic colours of ERGs, (^775 - K s ) > 4.0 (Vega). Based on 
the best-fit optical to near-IR SEDs, the large majority of 
the ERGs have been formed in short timescales (instan- 
taneously or with characteristic times r < lGyr). For 
many of them, the best-fit age is larger than the SED- 
derived time-formation scale, i.e. age > 1 Gyr. However, 
even when these characteristics would correspond to sys- 
tems undergoing passive evolution, the SED of most of 
these ERGs can only be properly fitted with consider- 
able amounts of dust extinction. We find that the me- 
dian of the colour excess in the 24 /im ERGs at z > 2 is 
E[B — V) « 0.30, a value too high to consider that these 
galaxies could be passive systems even from the optical 
to near-IR point of view. 

It should be emphasized that the SFRs computed here 
assume that the mid-IR sources are purely star-forming. 
Another possibility is that at least part of the mid-IR 
emission of these sources is due to the presence of an 
AGN, in which case the real SFR would be lower than 
our estimated values. We note, however, that even when 
e.g. only one half of the mid-IR luminosities of the z ~ 
2 — 3 galaxies analyzed here were due to star formation, 
the derived SFR would still be sufficiently high as to 
construct stellar masses of the order > 10 10 M Q in a 
burst timescale. The properties of these high-redshift 
massive mid-IR galaxies resemble those of submillimetre 
galaxies, whose redshift distribution has a median of z = 
2.2 - 2.4 ijChapman et al. 2 003: Ch apman et al. 2005[) . 
Thus, it is expected that the 24 /im galaxy population at 
the depth of the GTO images contains at least a fraction 
of the submillimetre galaxy population (cf. Egami et al. 
2004; Ivison et al. 2004). 

At lower redshifts z < 1, the completeness limits of the 
samples allow us to determine that star-formation spans 
a wide range of timescales, in agreement with results ob- 
tained from the analysis of ISO sources (Franceschini 
et al. 2003). Moreover, we find that these timescales 
for star-formation activity depend on galaxy mass, with 
more massive star-forming galaxies probably experienc- 
ing more prolonged star- formation histories. In particu- 
lar, for the high-mass end M > 1.5 x 10 n Af o at z < 1, 
we see that the derived timescales are larger than the 
age of the Universe at the corresponding redshifts. This 
fact indicates that the SFR of these galaxies was not 
constant over time, but larger in the past. This plausible 
differential star-formation history for star-forming galax- 
ies of different assembled stellar masses at z < 1.0 is in 
agreement with recent results reported in the literature 
(Heavens et al. 2004). 

Although the derived timescales Mj SFR for massive 
(M > 1.5 x 10 11 M Q ) galaxies at z < 1 indicate pro- 
longed star-formation histories, different possibilities ex- 
ist for the way in which star formation proceeded through 
cosmic time. Star formation could have proceeded quies- 
cently, i.e. with relatively small SFR sustained through 
long (> 0.1 Gyr) periods of time. Alternatively, as we 
discussed before, star formation might have been non- 
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continuous, with the production of multiple burst-like 
episodes. Both ways -or also some combination of the 
two- could have lead to the construction of the stellar 
mass contained in these massive galaxies. However, the 
number densities of both all and S v (24 /xm) > 83 fi Jy 
massive galaxies at z ~ 3 is only ~ 20% of the cor- 
responding number densities at z < 1 (Caputi et al. 
2005b, c). Thus, the major episodes of star- formation 
activity produced at very high redshifts could be respon- 
sible for the construction of stellar mass in at most a 
minor fraction of the massive 24 /im galaxies observed at 
z < 1. Also, it is likely that most of the massive galax- 
ies undergoing high star-formation activity at z ~ 2 — 3 
do not have significant amounts of remaining gas to con- 
dense into stars by redshift z ~ 1. Thus, massive 24 /im 
galaxies at z ~ 2 — 3 are probably not the progenitors of 
the massive 24 /im galaxies at z ~ 1. The star- formation 
history of massive IR galaxies at low redshifts should 
have mainly proceeded by some combination of quiescent 
activity and relatively modest burst-like episodes. 

It is interesting to analyze the evolution with redshift 
of the efficiency of the burst-like star-formation mode to 
construct a significant fraction of a galaxy stellar mass. 
As we said above, at redshifts z ~ 2 — 3, the typi- 
cal starburst lifetime of 0.01 Gyr to 0.1 Gyr is suffi- 
cient to form a substantial fraction of the stellar mass 
of some massive galaxies ( M > 1.5 x lO n M ). At 
redshifts z ~ 1 — 2, only low to intermediate-mass galax- 
ies (M < 7.0 x 10 10 Af Q ) could form in a burst lifetime. 
Within the limits of our sample and the volume surveyed, 
no galaxy appears to be able to have built up in a burst 
timcscale at z < 1. This does not mean that the burst- 
like mode cannot proceed in different mass galaxies at 
these low redshifts, but each individual burst would only 
be sufficient to construct an additional minor amount of 
the galaxy assembled stellar mass (i.e. the specific star 
formation rate SFR/M would be low). Thus, our results 
suggest that the potential importance of the burst-like 
mode of formation shifts from high to low-mass galaxies 
with decreasing redshift (note, however, that we cannot 
completely probe star-formation in galaxies with stellar 
mass M < lO n M at z ~ 2 - 3). In fact, each individ- 
ual burst-like episode can play only a secondary role in 
galaxy building at z < 1 (see also Papovich et al. 2005b) . 

6. DISCUSSION AND SUMMARY 

In this work, we presented _ftT s -band identifications of 
the MIPS 24/zm galaxies in the GOODS/CDFS at the 
depth of the Spitzer/ GTO surveys (with ~ 94% of iden- 
tification completeness for sources with flux S v (24 /im) > 
83/iJy). A K s < 21.5 (Vega) galaxy survey is sufficient 
for this purpose. The remaining non-identified galax- 
ies are mostly faint mid and near-IR sources. A minor 
fraction (7%) of the 747 identified objects correspond to 
X-ray-selected active galaxies. Active galaxies are among 
the brightest 24/im sources at high redshifts (cf. Houck 
et al. 2005). 

MIPS 24/tm galaxies in the GOODS/CDFS span the 
redshift range z = — 4 and 28% of the population lies 
at redshifts z > 1.5 at the depth of the GTO data. This 
confirms the predictions of a substantial population of 
mid-IR-selected galaxies at high redshifts. We deter- 
mined the existence of a significant bump in the redshift 
distribution of 24/im galaxies at z ~ 1.9, in agreement 



with predictions made by Lagache et al. (2004) using 
the Lagache et al. (2003) model. This secondary peak 
in the redshift distribution of 24/zm galaxies indicates 
the existence of PAH spectral features in z ~ 1.9 star- 
forming galaxies. We also observed the predicted depres- 
sion in the redshift distribution at z <~ 1.5, although with 
marginal significance, which could be driven by silicate 
absorption. However, the overall observed redshift dis- 
tribution differs significantly from the distribution pre- 
dicted by the Lagache et al.(2003) and other models, in- 
dicating that revisions are required to understand the 
high-rcdshift infrared-luminous galaxy population. 

The mid-IR galaxy population is mainly composed of 
normal star-forming galaxies at redshifts 0.4 < z < 0.8, 
while LIRGs start to be dominant at 0.8 < z < 1.2. We 
derived estimated assembled stellar masses of 10 9 M Q 
to 10 12 M© for the 24 /im galaxies at z < 1.2. Mas- 
sive star-forming galaxies appear to dominate the mid- 
IR output progressively at higher redshifts. A consider- 
able fraction of the stellar mass density of the Universe 
is contained in bright 24 /im sources: galaxies with flux 
S , „(24 /im) > 83 /iJy contain ~ 35% of the stellar mass 
density at redshifts z = 0.5 — 1.5 and ~ 40% of the stel- 
lar mass density at z ~ 1.5 — 3.0. This corresponds to ~ 
30%, 15%, 10% and a minimum of 5% of the local value 
at redshifts z w 0.75, 1.25, 1.75 and 2.5, respectively (Ca- 
puti et al. 2005c). Thus, we conclude that bright 24 /im 
galaxies have a significant role in the history of stellar 
mass assembly. 

We find that, at redshifts z < 1.2, the IR- 
derived timescales for star-formation activity increase 
with the already-assembled galaxy mass. Mid-IR 
galaxies with intermediate to high assembled stellar 
masses (M > 10 10 Mq) appear to be suffering rather 
prolonged star-formation histories, while lower-mass 
(M < 10 10 M Q ) galaxies, on the contrary, are charac- 
terized by star- formation activity on shorter timescales. 
At higher redshifts z > 2, massive galaxies are dominant 
and all the derived timescales for the S v (24 /im) > 83 /iJy 
sources are < 1 Gyr. Star-formation history at high red- 
shifts is likely to have proceeded in a series of burst-like 
episodes and, in a substantial fraction of massive galax- 
ies, a single burst lifetime (0.01 Gyr - 0.1 Gyr) is sufficient 
to assemble a stellar mass of w 10 10 M & - 10 11 M . We 
conclude, then, that the burst-like mode was probably 
an efficient way to construct the bulk of the stellar mass 
of massive galaxies at high redshifts. 

At z ~ 1 — 2, on the contrary, the burst-like mode can 
be effective only for the formation of galaxies with mass 
M < 7.0 x 10 10 M Q , and individual bursts can play only a 
secondary role in galaxy-building by redshifts z < 1 . Re- 
cently, Bell et al. (2005) found that the majority of the 
24 /im galaxies present at redshift z ~ 0.7 in the CDFS 
have regular spiral morphologies rather than irregular or 
peculiar morphologies. In conjunction, Bell et al. mor- 
phological determinations and our results on the lack of 
galaxies built up in a burst at z < 1 suggest that mor- 
phology could be linked to the mode of star-formation ac- 
tivity (see also Papovich et al. 2005b). Galaxies with no 
on-going or mainly quiescent star-formation would show 
regular shapes, while galaxies suffering important burst- 
like episodes of star-formation activity would appear as 
irregular. To test this hypothesis, a detailed morpholog- 
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ical study of 24 /im galaxies as a function of stellar mass 
and redshift would be necessary. 

Results recently published in the literature suggest 
a 'down-sizing' history for the formation of the stellar 
mass, where star-formation proceeds from high to low- 
mass galaxies (Heavens et al. 2004; Juneau et al. 2005). 
Our conclusions support a similar down-sizing mecha- 
nism for the efficiency of the burst-like mode. This fact 
also suggests that massive star-forming galaxies at high 
rcdshifts are not the progenitors of the massive star- 
forming galaxies found to be at z < 1, which are char- 
acterized by much longer star-formation histories. The 
properties of 24 /im galaxies at z > 2 suggest, on the 
contrary, a likely connection with submillimctre galax- 
ies. Current interpretations of submillimctre galaxies in- 
dicate that this population could contain the progenitors 
of the most massive E/SO galaxies present in the local 
Universe (e.g. Stevens et al. 2003). The high SFR char- 
acterizing mid-IR-selected galaxies at 2 < z < 3 suggests 
a similar conclusion for the high-redshift 24 /im galaxy 
population. However, some recent works suggest that 
a non-negligible fraction of present-day massive galaxies 
would be in place at even higher z > 3 redshifts, indicat- 
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Fig. 2. — 24 fim fluxes versus K s magnitudes for the MIPS 24 fim sources in the GOODS/CDFS associated to a K s < 21.5 counterpart. 
The circles, squares and asterisks correspond to normal galaxies, active galaxies and galactic stars, respectively. The dotted line indicates 
the limit of ~ 80% completeness of the 24 (im catalog, above which the if s -band identification is ~ 94%. 
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Fig. 3. — a) The redshift distribution of the MIPS 24 /im galaxies with _ft" s < 21.5 counterparts in the GOODS/CDFS: all sources and 
sources with S v > 83 fijy (empty and shaded histograms, respectively), b) The normalized redshift distributions of all and Sv > 83 fijy 
24 /im galaxies (black dashed and solid lines, respectively), compared to the normalized redshift distribution of the total K s < 21.5 sample 
in the same field (grey solid line). 




Fig. 4. — The 95% confidence limits on the normalized redshift distributions of the MIPS S„(24 fim) > 83 ^tjy galaxies (black lines) and 
all the K s < 21.5 galaxies (grey lines). The bump in the redshift distribution of 24 (im galaxies at z ~ 1.9 indicates the presence of PAH 
emission at these high redshifts. The marginally significant depression at redshift z ~ 1.5 might be due to silicate absorption. 
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Fig. 5. — 24 /im flux versus redshift for the MIPS 24 fim galaxies in the GOODS/CDFS associated to a K s < 21.5 counterpart. The 
symbols are the same as in figure [2] Filled symbols indicate spectroscopic redshifts. 



TABLE 1 

Completeness limits for the A" s -band 

IDENTIFICATIONS OF THE Spitzer MIPS 24 (im 
SOURCES IN THE GOODS/CDFS. 



24 fim flux 
S„ My] 


i<" s -band identification 
completeness 


> 250 


100% 


> 200 


99% 


> 150 


96% 


> 100 


95% 


> 83 


94% 
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Fig. 6. — The comparison of the observed rcdshift distribution of MIPS 24 (im galaxies in the GOODS/CDFS (solid histogram) 
with the redshift distribution predicted by Lagachc ct al.(2004) (dashed histrogram). Both distributions correspond to sources with flux 
S V (2A fira) > 83 /ijy in an area of 131 arcmin 2 . 
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Fig. 7. — The estimated assembled stellar masses for the MIPS 24 /im galaxies with K a < 21.5 counterparts in the GOODS/CDFS. The 
symbols with a cross correspond to galaxies with flux S V (2A (im) > 83 £tJy. The dashed line indicates the estimated mass completeness limits 
imposed by the K s < 21.5 survey. However, we note that these limits are basically irrelevant for the galaxies with flux Su (24 /im) > 83 /iJy, 
as the i-sT s -band identification completeness for these galaxies is ~94%. 
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Fig. 8. — Bolometric infrared luminosities versus assembled stellar masses for the MIPS 24 /im galaxies with K s < 21.5 counterparts 
in the GOODS/CDFS, in different redshift bins: a) 0.4 < z < 0.8; b) 0.8 < z < 1.2; c) 1.2 < z < 2.0; d) 2.0 < z < 3.0. The dotted lines 
delimit the region of completeness at the mean redshift of each bin, taking into account the S„(24 fim) = 83 fi Jy limit and the fact that our 
ii" s -band identifications arc almost complete above that flux (i.e. the mass completeness limits imposed by the K a = 21.5 cut are basically 
irrelevant above that limit). 
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Fig. 9. — The evolution of the star-formation rate of galaxies with different assembled stellar mass: a) M < 1.5 X 10 10 Mq; b) 
1.5 X 10 10 M Q < M < 7.0 X 1O 1O M ; c) 7.0 X 10 10 M Q < M < 1.5 X 10 n M Q ; d) M > 1.5 X 10 n M Q . In all the panels, the dashed line 
delimits the region of completeness of the observed star-formation rates as imposed by the flux limit S v (24fim) = 83/iJy. 
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Fig. 10. — The ratio between the assembled stellar mass and the instantaneous star-formation rate of MIPS 24 (im galaxies with K s < 21.5 
counterparts in the GOODS/CDFS, versus redshift. Different symbols indicate different assembled stellar mass ranges: M < 1.5 X 10 10 Mq 
(circles), 1.5 X 1O 1O M < M < 7.0 X 10 10 M Q (squares), 7.0 X 10 10 M Q < M < 1.5 X lO n M (triangles) and M > 1.5 X 10 n M Q (asterisks). 
The solid line shows the age of the Universe as a function of redshift, while the dashed line indicates an upper limit to the typical lifetime 
of a starburst (0.1 Gyr). 



